Finite average lengths in critical loop models 



Jesper Lykke Jacobsen and Jean Vannimenus 

Laboratoire de Physique Statistique^, 
Ecole Normale Superieure, 
24 rue Lhomond, F -15231 Paris Cedex 05, France. 

ABSTRACT 

A relation between the average length of loops and their free energy is obtained 
for a variety of 0(n)-type models on two-dimensional lattices, by extending to finite 
temperatures a calculation due to Kast. We show that the (number) averaged loop 
length L stays finite for all non-zero fugacities n, and in particular it does not 
diverge upon entering the critical regime (n — > 2 + ). Fully packed loop (FPL) 
models with n = 2 seem to obey the simple relation L = 3L m i n , where L m ; n is 
the smallest loop length allowed by the underlying lattice. We demonstrate this 
analytically for the FPL model on the honeycomb lattice and for the 4-state Potts 
model on the square lattice, and based on numerical estimates obtained from a 
transfer matrix method we conjecture that this is also true for the two- flavour 
FPL model on the square lattice. We present in addition numerical results for the 
average loop length on the three critical branches (compact, dense and dilute) of 
the O(n) model on the honeycomb lattice, and discuss the limit n — * 0. Contact is 
made with the predictions for the distribution of loop lengths obtained by conformal 
invariance methods. 

PACS numbers: 05.50.+q, 64.60.-i, 75.10.Hk 
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1 Introduction 



Loop models appear in various contexts in statistical physics, not only as models of 
various phases of polymers [0 but also in the high-temperature expansion of the 0(n) 
spin model [Q, |3|, in surface growth models |3J], in efficient cluster-flipping Monte Carlo 
methods ||, and in the study of quantum spin systems ||. They have deep connections 
with other well-studied statistical models, such as the eight-vertex, colouring and Potts 
models |?|, |[, and they can be viewed as toy models for the study of string theories. One 
of their main interests is that they provide examples of systems involving extended degrees 
of freedom, rather than local ones as in spin systems, for which many exact results have 
been obtained in d = 2 using Bethe ansatz J7| , Coulomb gas || and conformal invariance 
methods [TO]. 



The basic O(n) loop model is defined through the partition function 

Z(n,t)=Y,n P t v , (1.1) 
c 

where n is the loop fugacity, which controls the number P of different loops ("paths") 
in the configuration C, and t is a temperature-like variable which controls the number 
V of empty sites ("voids") not covered by any loop, the total number of sites on the 
lattice being iV.Q For non-crossing loops on the honeycomb and square lattices a phase 
transition occurs at n = n c = 2, for all temperatures below a lattice-dependent t c || [lTJ. 
That transition is very weak and is reminiscent of the Kosterlitz-Thouless transition for 
the XY spin model, as the free energy only has an essential singularity. The geometrical 
interpretation is simple: For n > n c the configurations contributing to Z only contain 
finite loops in the thermodynamic limit, whilst for n < n c there appear "infinite" loops 
whose extent diverges with the lattice size. This is analogous to what happens in lattice 
percolation models, where an infinite cluster appears when the site (or bond) occupation 
probability is above a threshold value p c ||12|| , and indeed the perimeters of percolation 
clusters can be viewed as special cases of loops [JTO], |T5 |. 



2 This is the standard definition of the O(n) model on the honeycomb latttice ||. On lattices with 
higher coordination number a given site may be visited by more than one loop, necessitating the intro- 



duction of further weights in Eq. (1.1) [111. We shall come back to this point later on. 
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Explicit formulae for the free energy per site, defined as 



F(n,t)= lim i-logZ(n,t), (1.2) 

N—+00 jV 



can be obtained |14], |15|, [16| along certain critical lines in the (n, t) plane || for the 



honeycomb lattice, but the distribution of loop lengths cannot be inferred from the 
partition function, and information about it is obtained numerically |§ or through scaling 
and conformal invariance methods |17| (see the lectures by Cardy |10[ for an introduction 
to the geometric properties of loop models). We show in the first part of the present work 



that a simple relation, first derived by Kast at t — in the honeycomb O(n) model [[18 
can be generalised to finite temperatures and to a variety of other loop models. It gives 
the average loop length L from the knowledge of the free energy and its derivatives with 
respect to n and t only. 

In the second part we use the exactly known results for F(n, t = 0) to calculate ana- 
lytically that average length at the fugacity where a transition occurs, n c = 2. We present 
next very accurate numerical calculations of the average loop length, using a transfer ma- 
trix method and our relation, at finite temperatures and for the zero-temperature 0(n) 
model on the square lattice. For both the t = (fully packed) model on the honeycomb 
lattice and for the 4-state Potts model on the square lattice we find the exact result 

L = 3L m ; n , (1-3) 

where L m ; n is the minimal loop length allowed by the given lattice. Contrary to what was 



suggested in Ref. |18[ on the basis of numerical results, the average length thus remains 
finite for n = n c , and indeed for all non-zero n. This is a priori counter- intuitive, as 
one would naively expect the appearance of infinite loops below n c to be reflected in a 
diverging average length, but we show in the last part of the paper that this result is not 
in contradiction with the predictions of conformal invariance. 

2 Average loop length and the free energy 

The average length of non-crossing loops in a given configuration C on the honeycomb 
lattice is given by 

N — V 

(L{C)) = — pA (2-1) 
2 



since a given site belongs to at most one loop, and its ensemble average with respect to 
Eq. Q is 

I = ivt!„V. (2 . 2) 

Z 1 ^ P v ; 

Introducing the partial derivatives d n and dt with respect to n and £ one has: 

d n (ZL) = l^(N-V)n p t v , (2.3) 
n c 

d t Z = j^y/f. (2.4) 
Eliminating the average number of voids between these two relations gives 

d n (ZL) = ±-[NZ -td t Z], (2.5) 

or, in terms of the free energy, 

^ + Ld n F=^[l-td t F]. (2.6) 

Now, if d n L remains finite the first term on the left-hand side of Eq. ( |2.6j ) becomes 
negligible in the thermodynamic limit, and we obtain a simple relation between L and 
the partial derivatives of F(n,t), 



which generalizes Kast's relation [18| to non-zero temperature. 



Several remarks are in order here: 

• The validity of relation ( |2.7| ) relies on the finiteness of d n L, which in turn, as is 
seen from Eq. (j2.7|) , holds as long as d n d t F and d\F remain finite, and n and d n F 
do not vanish. As we shall see below, these conditions can be shown explicitly to 
be fulfilled for fully packed (t = 0) loops on the honeycomb lattice, for all n > 0. 

• The knowledge of F(n, t) along a critical line in the (n, t) plane is not sufficient to 
obtain the average length L, as it provides only one relation between d n F and %F, 
through the total derivative d n F, and in general it is not possible to express L as 
a function of d„F alone. 
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One can also write relations for P, the ensemble average number of loops [ H) : 



P = ^T,P(C)n P t v . 



(2.8) 



One has from Eqs. ( |1.1|) and ( |Q|) 



77 

P = -d n Z = Nnd n F. 



Taking now the derivative of Eq. ( |2.9| ) with respect to n, one finds 

2" 



(2.9) 



d n P 



^ d„Z 



Z 



n 



mz o n z 



z 



z 



-(p^-p 2 

n v 



(2.10) 



which yields an explicit expression for the fluctuations in the loop number: 

1 



N 



p2 _ =nd n F + n 2 d'F. 



(2.11) 



Relations (|2.9j ) and (|2.11 ) show that the average number of loops is extensive (ex- 
cept for n = 0) and that its fluctuations become negligible in the thermodynamic 
limit, as long as d n F and d 2 F remain finite. It is also interesting to note that 
for a perfect gas at fixed chemical potential the fluctuations in the number M of 
particles in a given volume are given by 



M 2 — M = M, 



(2.12) 



so the second term on the right-hand side of Eq. ( p.ll|) may be viewed as a correction 
with respect to the fluctuations in a "perfect gas" of loops, induced by the non- 
overlapping condition. 

Before leaving this Section we briefly comment on the appropriate form of the re- 
lation ( |2.7| ) for lattices with higher coordination number. As an example, consider the 
O(n) model on the square lattice, for which the partition function has a slightly more 
complicated form than the one given in Eq. ( [[ . 1|) , since each vertex can be visited by the 
loops in several ways that are unrelated by rotational symmetry. Following Ref. [[□]] we 
define 



c 



(2.13) 
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where N t , N u , N v and N w are the number of vertices visited by respectively zero, one 
turning, one straight, and two mutually avoiding loop segments. The relation ( |2.1|) must 
be replaced by 

WO) = *^±* (2 , 4) 

and under assumptions on the higher derivatives identical to those given above the gen- 
eralisation of Eq. ( |2.7|) finally reads 

Ld n F = -(l-td t F + wd w F). (2.15) 
n 

3 Zero-temperature case: Fully packed loops 

At t = configurations containing voids do not contribute to Z, and on the honeycomb 
lattice every site must belong to one and only one loop. The system of loops is then 
called fully packed or "compact" , to distinguish it from the "dense" phase 0, [T7]] with a 
finite density of voids found at non-zero t. A necessary condition for the t = line to 
be stable (in the RG sense) is that Z(n,t) is invariant under the replacement of t by —t 
||19||f| This symmetry criterion is fulfilled because the honeycomb lattice is bipartite and 
hence only allows for loops of even length. For the model at hand the condition is also 
sufficient (see Ref. |2(J for an example where this is not the case), and as a consequence, 



the t = line is critical for n < 2 and fully packed loops are in a different universality 



class than dense loops |19|, ^T], [22|, |23], |4|. The relation (|2.7p reduces for t = to the 
form obtained by Kast |18| : 

Lin) = }, r , (3.1) 
v ; ni^(n)' v ; 

where F (n) = F(n,t = 0), and it is particularly interesting as explicit expressions are 
known for the free energy in that case, thanks to the work of Baxter [J25|J , Reshetikhin 
26| , Batchelor et al. | 21| , and Kast himself |Tj|. One has for the free energy per site of 



the honeycomb lattice 
• for n < 2: 



„ . . 1 r°° , sinh 2 Xx sinh(7r — X)x 



where A > and n = 2 cos A. 
3 Negative values of t correspond to an antiferromagnetic interaction in the spin language. 
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for n > 2: 



F o( n ) = 2 lo § 



« 1/3 n,7T 



-6m+2\2 



1 (1 -g- 6m+4 )(l -g" 6m ) 



(3.3) 



where g = e 7 and n = 2 cosh 7. 



Analytically it is known that at the critical point n c = 2, F (n) takes the value |2f| |2l| P| 



*b(2) 



du sinh u 



m sinh 3u 



2 00 

z m=l 



(3m - l) 2 



3r 3 (i/3) 

47T 2 



0.18956' 



(3.4) 



The analytic calculation can be pushed further, noting that for n < 2 Eq. ( |3.2| ) can be 
reexpressed as 



F (A) - F (A = 0) = / 

Jo 



00 sinh u 
du 



u sinh 3m 



1 - 



1 



tanh(7ra/A) 



00 dv sinh 3 (At>/7r) 



v sinh(3Af/7r) 



tanhw 



(3.5) 



This expression can then be expanded in powers of A. To lowest order this reads: 



F (A) - F (A = 0) 



A 5 



dvv 



2A 2 r°° , v 
dv 



3tt 2 Jo sinh v 3n 2 Jo e 2v — 1 

In terms of the loop fugacity n, this yields 



dF n 



dn 



dFndA 



dA dn 



1 

36' 



X 2 



(3.7) 



so we obtain an unexpectedly simple result for the average length of the fully packed 
loops at the critical point 

L (n = 2) = 18. (3.8) 



A similar analysis can be performed for n > 2, rewriting for example Eq. ( |3.3j ) as 

7 1 



F (n) - F (2) = ± - ^ £ /„ 

6 2 m=l 



with 



fm = log 



g(6m7)g([6m — 4)7) 
2 2 ([6m-2] 7 ) 



(3.9) 



(3.10) 



Note that there is a misprint in Eq. (7) of Ref. fell 
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and g(x) = ~(1 — e x ). The terms in Eq. (|3.9| ) have been grouped so as to make it more 
convenient to use the Taylor-McLaurin summation formula 

00 roo 11 1 

Ef™ = J o di//(y)--/(0)--/'(0) + — r(0) + .-.. (3.11) 

The evaluation of Eq. (3^9) for small 7 then yields for dF /dn\ n ^ 2+ the same value as 
Eq. 0. 

The calculation of d 2 F /dn 2 proceeds along similar lines, and the end result reads 

d 2 F 1 



dn 2 



1080 



for n 



(3.12) 



the same value being obtained on both sides of the transition (more generally, one expects 
all derivatives of F (n) to be continuous for n — 2). As A = (n = 2) is the only point 
where Eq. ( |3.2|) may have a singularity, we conclude that Kast's relation (|2.7|) remains 
valid in the whole region n < 2, which was not a priori obvious from his derivation. 

As a by-product one also obtains from Eqs. ( |2.9|) and ( 2.11] ) the average number of 
fully packed loops for n = 2: 



n=2 



N 
18 



and its fluctuations: 



Po 2 ~ Po 



n=2 



(3.13) 



(3.14) 



.18 270, 

Another interesting case is the limit of n —>■ 0. This is the compact polymer (Hamilto- 
nian cycle) limit, when a single loop covers all the lattice sites (with adequate boundary 
conditions). The derivative of Eq. (|3.2|) can now be evaluated directly: 



dF 



dA 



Thus 



X=n/2 

dF 



4 

7T JO 



00 , 1 - cosh 2u 
du 



dn 



l + 2cosh2w) 2 vr 3V3' 
dF dA 



n=0 



1 



1 



A=tt/2 



3V3 2tt : 



dA dn 

and invoking once again Eq. ( |3.1| ) the average loop length L (n) diverges as 

30.0344 • • • 



(3.15) 



(3.16) 



Loin) 



n 



0(1). 



(3.17) 



One can also show that the second derivative F"(0) = —0.00684676 • • • is finite. 
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4 Numerical results on the honeycomb lattice 



4.1 Zero temperature 

We have checked numerically the analytical expressions given in the literature for the free 
energy of the O(n) model at t = 0, using a previously developed transfer matrix method 
T9fl to compute the partition function for strips of finite width, up to W mSuX = 15. The 



results can be further analysed by exploiting that the dominant finite-size corrections are 
known from conformal invariance arguments E7I, |2£ 



Fo(WO = Fo(oo)-^ + ..., (4.1) 
where c(n) is the central charge of the model. The latter has been inferred from the 



Bethe ansatz solution of the model ]2l| 



6e 2 

c(n) = 2 with n = 2 cos(7re). (4.2) 

1 — e 

Carrying out the numerical differentiation on Fq(oo) obtained from Eq. (|4.1|) , rather than 
on the Fq(W) which are related to the leading eigenvalue of the transfer matrix spectra, 
we were able to obtain sufficient accuracy to calculate the derivative dFo/dn, and hence 



Lo(n), using formula (|2.7|) . 

In Fig. |l] we plot Lo(n) in units of L m m = 6, which is the shortest possible loop 
length on the honeycomb lattice. To emphasize the fact that the average loop length 
does not exhibit any singularity at n c = 2 these data are shown without the central 
charge correction. We point out that it would not have been possible to produce this 
plot by a direct numerical differentiation of the analytic results, Eqs. ( |3.2| ) and ( |3.3| ). 
as these expressions in their present form converge extremely slowly close to the critical 



fugacity n c = 2. Actually we suspect that it was these difficulties that led Kast [18| to 
the erroneous conjecture that Lin — > 2 + ) diverges. 

Turning now on the correction (|4.1|) we find the finite-size values of Fq(u) given in 
Table |I[ Results are only shown for W a multiple of three, since otherwise the equivalent 



interfacial representation |22| will be subject to height defects, leading to the introduction 



of non-trivial twist-like operators in the continuum theory [p3[| . An extrapolation to the 



thermodynamic limit W — > oo can be performed by fitting the residual size dependence 
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1 2 3 4 5 

r 



Figure 1: Average loop lengths, measured in units of six, for the FPL model on the 
honeycomb lattice. The finite-size estimates converge rapidly towards the exact value 
L(n = 2) = 18. 



w 


F>(n = 0) 


6 


0.03192052 


9 


0.03304791 


12 


0.03321945 


15 


0.03326456 


oo 


0.033296 (3) 


Eq. (|3.16|) 


0.0332951- •• 



Table 1: Finite-size estimates of F^(n = 0), corrected for the known value of c, as a 
function of strip width W. The extrapolated value is found by fitting the residual size- 
dependence to a 1/W A term. 
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to a 1/W 4 term [23]. The result is 



30.034(3), 



(4.3) 



F^n = 0) 

in excellent agreement with the analytical result ( 3.17|) . In the same manner we find for 



n = n c that 



35.96(3), 



(4.4) 



K(n = 2) 
confirming Eq. fl3.7| ). 

We have also evaluated numerically Fq(ji = 0) = —0.00686(2), whereas an attempt to 
compute Fq{7i = 2) did not lead to a sufficiently accurate result to permit a comparison 
with Eq. ( |3.12| ). The reason for this is a combination of the logarithmic corrections to 
the scaling form 1|) |29| and the singularity arising in Eq. (|4.2|) when taking the second 
derivative with respect to n at n = n c . 



4.2 Non-zero temperature: Dense and dilute loops 

As discussed earlier, explicit expressions for F(n,t) are known only along two other 
special critical lines in the (n, t) plane, for which Bethe ansatz solutions were given in 
Ref. |16|]. These lines, or branches, are given by the relation |J 

t 2 = 2 ± V^n, (4.5) 

where the upper sign corresponds to the "dilute" branch and the lower sign to the "dense" 
branch. We have already remarked that this knowledge is not sufficient to obtain analyt- 
ically the average loop length from Eq. ( |2.7| ). However, an accurate numerical determi- 
nation of the partial derivatives d n F and d t Ft is possible, using the same transfer matrix 
approach as for t = so the variation of L along these lines can be studied in detail. The 
central charge values used in the extrapolation procedure are now 

6e 2 

c(n) = 1 with n = 2cos(7re). (4-6) 

1 =F e 

The resulting values of L(n) are displayed in Fig. [2|. Only the data for strip width 
W = 9 are displayed, since employing the correction ( |4.6| ) renders the residual finite-size 
effects indistinguishable on the figure. Extrapolating to W = oo we obtain for n = n c = 2 

L(nc,t c = V2) = 19.417(1), (4.7) 
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Figure 2: Average lengths of dense and dilute loops on the honeycomb lattice, measured 
in units of six. 

which is slightly larger than the corresponding value at t = 0. The same is true for the 
asymptotic n — > behaviour on the dense branch, for which the numeric result is 

4n S c(n) = ^ + 0(l). (4.8) 
n 

At first sight this finding may seem odd, since the geometrical scaling dimension 

x > = ¥ - e) - 2{&7i < 4 ' 9) 

coincides for the dense and the compact branch R BT[, implying that the asymptotic 



distributions of large loops are identical in the two models [17]. However the dense phase 
has fewer short loops than the compact one, since it is always energetically favourable to 
trade a length-six loop for some empty space: 

t 6 (2 - V / 2 T7 ^) 3 



n n 



> 1 (identity at n = 1). (4.10) 



Since a length-six loop necessarily has a fixed shape of its perimeter there is no extensive 
entropy involved in the argument, and our simple energetical consideration suffices. 

Another remark pertains to the fact that L(n) is an increasing function of n on 
the dilute branch (and in particular it does not diverge as n — ► 0). This has a quite 
elementary explanation in terms of X2 to which we shall come back in Section |6|. 
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5 Loops on the square lattice 



Many results have also been obtained for systems of loops on the square lattice [0, || |TT 



23| , |24] , ftO[ |3T| , |32| , p3| , |34|, but now there are several possible definitions of the models, 
depending on whether one allows the non-overlapping loops to cross at the vertices or 
not, if there is only one kind of loops or two, and exactly how the loop segments are 
allowed to turn at the vertices. Following Refs. ]23|, p4j| , we consider a model with two 
kinds ("flavours") of loops, with respective fugacities n\ and n 2 , which both cover all N 
vertices of the square lattice and can intersect but not overlap. This can be considered 



as a particular zero-temperature limit [[34] of the O(n) model on the square lattice [|TT 



The partition function and the free energy are given by 

Z (n 1 ,n 2 ) = J2 n i ln 2 2 i 
c 

F (n 1 ,n 2 ) = Jim — logZ (rn, n 2 ), (5.2) 

where the summation is performed over all doubly compact configurations. The con- 
straint that both types of loops cover all N sites implies that for each configuration the 
average length of type-1 loops is just (Li(C)) = N/Pi, and its ensemble average is given 
by 

^=^E^< 1 < 2 - (5-3) 
Repeating the calculation of Section |^ one obtains 

d ni (Zj7 l ) = N^- (5.4) 

m 



and 



i^Lx + LA.Fo^— . (5.5) 
iV rii 



Now, if d ni Li remains finite, the first term in the left-hand side of Eq. ( |5.5| ) is negligible 
in the thermodynamic limit, and we obtain 

Ll nid ni Fv 

A similar relation holds for the average length L 2 of the second kind of loops, with n 2 
replacing m, and when ri\ = n 2 both reduce to Kast's relation ( [3.1|) . 
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One can also define the average length L of any loop, regardless of flavour. Following 
the above line of reasoning this is easily shown to fulfill the relation 

L = 2=EL. (5.7) 
U + L 2 K J 

For n% = n 2 the above statement is trivial. Taking however the limit L 2 — > oo we obtain 
the curious result that 

L = 2L[ for n 2 -> 0, (5.8) 

regardless of the value of ri\. 

Since the model ( |5.1|) has as yet defied a Bethe ansatz solution, no analytical expres- 



sions for ^0(^1,712) are available (apart from the conjecture i*o(0, 0) = |log2 ||23| ). We 
have therefore studied numerically the relation (|5.6| ) as a function of ni, for different 
values of n 2 — 0, 1 and 2. Transfer matrix computations were performed for strips of 
even width, up to W max = 12, exploiting that |2"3[ 

6e? 



c(ni,n 2 ) 



ei 



6e| 
l-e 2 ' 



(5.9) 



Li(ni = 2, n 2 ) 



^4 

L x {ni -> 0, n 2 ) = — + C(l) with 

71^ 



Extrapolating as usual to = 00 we find for the average loop lengths at rt\ = 2: 

15.74(5) for n 2 = 

13.95(3) for n 2 = 1 (5.10) 
11.99(2) for n 2 = 2 

and for its divergence when rti — » 0: 

A = 21.97(6) for n 2 = 
A = 20.35(3) for n 2 = 1 ■ (5.11) 
A = 18.96(3) for n 2 = 2 

A striking fact is that for the symmetric case rti = n 2 = n c = 2, which corresponds to 
the four-colouring model introduced in Ref. 0, the numerical value of the average length 
of either kind of loops is 

T c /L min = 2.997(4), (5.12) 

where we have normalised with respect to the length of the shortest possible loop on the 
square lattice, L min = 4. The ratio (|5.12|) is extremely close to 3, which by Eq. (|3.8| ) is 
the exact result for the zero-temperature O(n) model on the honeycomb lattice. 
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In order to see if this could be more than a coincidence, we have looked at another 
model on the square lattice, describing the g-state Potts model at its selfdual point [0. 
We briefly recall how the Potts model is transformed into a loop model M : In a first step 



the spin model (with coupling constant K) is turned into a random cluster model [|35 
Each state is now a bond percolation graph in which each connected cluster is weighed 
by q and each coloured bond carries a factor of (e^ — 1). Next, the clusters are traded 
for loops on the medial graph (the lattice formed by the mid-points of the bonds on the 
original lattice). By the Euler relation each closed loop is weighed by ^/q, and the bond 
weights are (e K — l)/^/? [0]. At the selfdual point the latter equals one by duality, and 
the exact correspondence is 

Z Potts = q N/2 T,V^ P = 1 N/2 Zloo P , (5.13) 
c 

where N is the number of sites of the original lattice. 

Defining F as the free energy per medial site with respect to Z^ OOV) one has |7j 



for n < 2: 



. 1 r°° , sinh Ax sinh(7r — X)x 

Fin) = - dx v - - — , (5.14) 

2 J-oo x smh 7rx cosh Ax 



where A > and n = ^/q = 2 cos A. 
for n > 2: 



A °° 6~ m ^ 

F{n) = tanhmA, (5.15) 

2 m=i m 



where n = ^fq = 2 cosh A. 

The value at the critical point n c = 2 is also known: 

du 



poo Hi), 

F(n = 2)= — tanhwe"" = 0.78318- ••, (5.16) 
Jo u 



and as in Section [3] we can evaluate the derivative dF/dn| n ^ 2 _ by examining the differ- 
ence 

roc fa 



F(X) - F(X = 0) = - r — sinh(At;/7r) tanh(Aw/7r) 

Jo v 



1 -i 



tanh v 



(5.17) 



At the bottom line we obtain 

L(n = 2) = 12, (5.18) 
14 



or L c / L min = 3, the same result being found also for n — > 2 + . For n — > the average 
loop length diverges as 

L(n) = — + 0(1). (5.19) 

n 

Both these results are in excellent agreement with numerical estimates produced by the 
transfer matrices described in Ref. [36[ . 



Our computations for the Potts model suggest that the result ( |1.3|) is not a mere 
coincidence, and we feel confident in conjecturing that L/L m i n = 3 is also exactly true 
for the four- colouring model, cf. Eq. ( |5.12| ). It would be quite surprising if this ratio were 
a new universal quantity, like for instance the ratio of the average loop area to the average 
square gyration radius, calculated by Cardy JTOf- Rather it is reminiscent of the "quasi- 
invariants" found in percolation theory: For a large variety of lattices, the percolation 
thresholds are observed to be little lattice-dependent, for a fixed dimensionality, once 
they are properly normalised taking into account the coordination number |37] . 



6 Discussion and contact with conformal invariance 
predictions 

The most striking aspect of the results presented above is that the average loop length 
remains finite, even inside the critical phase where the correlation length, which is roughly 
the size of the largest loop in a typical configuration, diverges. In fact, a little thought 
shows that the apparent paradox comes from the definition of "average length" : In our 
calculations all loops are equally weighted, irrespective of their length. This is called 
the "number average" in polymer theory, to distinguish it from the "weight average", 
where each polymer contributes with a weight proportional to its length. Depending 
on the quantity being measured, these and still other averages are relevant in order to 
compare theoretical predictions with experiments (see, e.g., the treatise by des Cloizeaux 
and Jannink for a detailed discussion). An analogous situation arises in percolation 



theory, when one defines the "average cluster size" |12j : If one means "the average size of 
the cluster on which an 'ant' lands at random", the weight average has to be taken, which 
is different from the average size obtained by simply giving every cluster the same weight. 
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In practice these various averages may be very different when the quantity considered has 
a power law distribution, as the result may be dominated in one case by the few largest 
polymers or clusters, and in another case by the numerous small ones. This is precisely 
what happens in the loop systems considered here: In the critical phase (n < 2) the 
(ensemble averaged) probability distribution of loop lengths (unweighted) is expected to 
be of the form 

II(L) ~ L~\ (6.1) 



where the exponent r is related to the geometrical (string) scaling dimension [17, || x 2 
by 

t-1 = -2— . (6.2) 
2 - x 2 

Conformal invariance arguments have permitted the exact evaluation of x 2 for all loop 



models discussed in this paper. These are P, O, El], 23, 24 



1 e 2 

x 2 = -(1 =F e) - 2 ^ ^ with n = 2cos(7re), (6.3) 

where the upper sign holds for any one of the compact or dense phases, and the lower 
sign for the dilute O(n) phase (see Ref. |RJ for a discussion of this kind of universality). 



In particular 3 > r > 2 for all < n < 2. This implies that the number averaged loop 
length 

E LU(L)\/[ E (6.4) 

^ = ^min / \^' = ^'min / 

converges, as we have shown explicitly. On the other hand, one can pick a point at 
random on the lattice and ask what is the average length of the loop to which it belongs. 
That weight averaged length is given by 

^max \ / ^max \ 

E L 2 U(L)\/[ £ LU(L) J (6.5) 

and it diverges like for large system sizes. In these equations L max is not the largest 
loop that can be placed on the iV-site lattice (a Hamiltonian cycle of length N), but 
rather loop size up to which the scaling law ( |6.1| ) is valid. On an i x i lattice this is given 
by L max ~ £ D{ , where |17| 

D { = 2 - x 2 (6.6) 



is the fractal dimension of the loop. Using Eq. ( |6.2j ) we then find 

L* ~ (£ 2 ) D{ - 1 ~ N Dl ~\ (6.7) 
16 



Figure 3: Detail of a typical configuration in the four-colouring model. The two loop 
flavours consist of the black and white lattice edges respectively. A loop passing through 
a randomly chosen point is likely to be very long, like the one of length 1076 shown in 
boldface. The presence of numerous short loops, however, explains why the (number) 
averaged loop length is finite in the thermodynamic limit: L = 12. 

For the four-colouring model Kondev and Henley M have confirmed this picture in 
detail using Monte Carlo simulations with non-local loop updates. In particular the 
power-law distribution of loop lengths ( |6.1| ) was confirmed very accurately, with the ex- 
pected value t = 7/3. However, these authors did not examine the finite-size dependence 
of L*. In order to verify the validity of Eq. (|6.7|) we have performed Monte Carlo simula- 
tions along the lines of Ref. || on several lattice sizes, up to £ max = 200. In each case the 
system was equilibrated by means of 10 5 loop flips, whereafter the lengths of a further 
10 6 loop updates were registered. (A detail of a typical loop configuration is shown in 
Fig. |.) We found that 

in good agreement with Eq. ( |6.7| ), since the fractal dimension is known to be exactly 
D { = 3/2 ||. Unfortunately, it is not possible to relate L* directly to the partition 
function and to obtain relations analogous to those holding for L. 

Let us end by commenting on the divergence of L as n — > 0. By naively integrating 
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Eq. ( |6.1|) with a lower cut-off L min and using scaling relations one finds 

9 T . 

L( n ) ~ (6.9) 
x 2 

Since from Eq. ( [4.9| ) x 2 — > as n — > 0, it is not surprising that L is found to diverge 
for compact and dense loops. On the other hand, for the dilute phase x 2 as given by 
Eq. fl6.3| ) does not vanish as n — > 0, in accordance with Fig. In fact one has the rough 
estimate 

(WO) 

L(n = 0) 4 v ; 

which is not too far off the numerical value of ~ 0.72. 
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